Mean field phase diagrams of imbalanced Fermi gases near a Feshbach resonance 
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We propose phase diagrams for an imbalanced (unequal number of atoms or Fermi surface in two 
pairing hyperfine states) gas of atomic fermions near a broad Feshbach resonance using mean field 
theory. Particularly, in the plane of interaction and polarization we determine the region for a phase 
separation phase composed of normal and superfluid components. We compare our prediction of 
phase boundaries with the recent measurement, and find a good qualitative agreement. 
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Two recent experimental studies of fermionic super- 
fluidity in strongly interacting atomic 6 Li gases with 
controlled population imbalance in two spin compo- 
nents have attracted intense interest from the physi- 
cists in wide communities 0, 0|. A very salient rea- 
son is the mysterious nature of the pairing mechanism 

SSESHSIiESIIIi- Since the Bardeen-Cooper- 
Schrieffer (BCS) pairing requires an equal number of 
atoms in each spin state, the presence of spin popula- 
tion imbalance leads to some exotic forms of pairing, 
such as the finite-momentum paired Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) state the breached pairing or 
Sarma superfluidity and phase separation jfj- How- 

ever, the true ground state of imbalanced fermionic su- 
perfluidity remains elusive and has been the subject of 
debate for decades. The two recent experimental obser- 
vations open up the intriguing possibilities for resolving 
this long-standing problem. As the population imbalance 
increases, the disappearance of superfluidity has been 
identified and phase separation of a unitary gas in 
trap has been observed Q. 

Motivated by the significant experimental develop- 
ment, in this paper we present a general mean field 
analysis of the ground state of homogeneous imbalanced 
atomic gases, focusing on the strongly interacting region 
near the broad Feshbach resonance, namely, the so-called 
crossover from BCS superfluidty to the Bose-Einstein 
condensation (BEC). Our goal is to map out the qualita- 
tive zero-temperature phase diagrams in the entire BCS- 
BEC crossover. A previous discussion of such phase di- 
agrams is based on a purely educated guess Q. Further 
analytic mean-field estimate is restricted to the narrow 
Feshbach resonance @, for which the most fascinating 
crossover region has been essentially ruled out, and thus 
is of less experimental relevance. 

In contrast to these prior theoretical studies, our anal- 
ysis is in close connection to the experiment and has more 
predictive powers. Our main results may be summarized 
as follows: (1) Aside from the ability to include the ex- 
otic phases mentioned earlier, our mean-field calculation 
predicts a new phase (the saddle point solution below), 



which becomes energetically favorable for a finite popula- 
tion imbalance. However, the new solution is inherently 
unstable towards phase separation, signifying an inho- 
mogeneous mixed phase. Around the crossover, in con- 
sistent with the experimental observations [l|, l2j we find 
that the phase separation phase becomes dominant in the 
phase diagram. (2) We construct the phase boundary of 
superfluid-to-normal transitions, and compare it with the 
measurement by Zwierlein et al. Q. The agreement is 
qualitatively good. 

We consider an imbalanced Fermi gas of 6 Li atoms 
across a broad Feshbach resonance, which is well de- 
scribed by using a single-channel model 
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Here the pseudospins a =|,| denote the two hyperfine 
states of 6 Li, and c^ a is the fermionic creation operator 
with the kinetic energy £ko- = £k — and e k = ft 2 k 2 /2m. 
The chemical potentials are different, i.e., jUfj = fi±Sfi, 
to account for the population imbalance 8n = n-j- — nj.. g 
is the bare interaction strength, and is expressed in terms 
of s-wave scattering length a via (Awfi 2 ! a / 'm) -1 = g + 
Ek(2e k )- 1 . 

In the mean field approximation we decouple the inter- 
action term by introducing an order parameter of Cooper 
pairs in momentum space A = — g^2 k (c q /2-ki c q/2+kT)i 
where the pairs may possess a nonzero center-of-mass mo- 
mentum q in case of spatially modulated states • As 
a result, the order parameter in real space aquires a one- 
wave oscillation form: A(x) = — g (c;(x)cf(x)) = Ae lq ' x . 
The value of q, together with A, are to be determined. 
The Hamiltonian can then be approximated by, 



H = ^£k<rc£ CT c k(r - [c|_ ki C| +kT + /i.C. 
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where in the second line we define a Nambu creation field 
operator: tp£ = (Cq /2+kT , c q / 2 -ki), and a z are the 2 x 
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Figure 1: (color online). Landscape of the thermodynamic 
potential at l/kpa = —1. The chemical potential is fixed to 
ji = 0.98942eF. The competing ground states are (i) a normal 
Fermi gas with A = 0, (ii) a fully paired BCS superfluid with 
A = Ao, q = 0, and Sn = 0, (iii) a breached pairing or Sarma 
superfluid with A < Ao, q = 0, and Sn / 0, (iv) a finite 
momentum paired FFLO superfluid with A < Ao, g / 0, and 
Sn / 0, and (v) a saddle point phase intervening between the 
local BCS and FFLO minima. 
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Figure 2: (color online). Phase diagram in the plane of inter- 
action and chemical potential difference. 



2 Pauli matrices, £k± = (£ q /2+kt ± £q/2-kj.)/2, and £q = 
Sk (? k + — Ck-) — A 2 /g. The above pairing Hamiltonian 
may be solved by the standard Bogoliubov transforma- 
tion, or more straightforwardly, by employing the Nambu 
propagator G(k, iu m ) = l/[(iu m - £ k -) - 6c+fz + Ao>] 
with quasiparticle energies Ek± = (^ + + A 2 ) 1 / 2 ± £k-. 
Here uj m = (2m + 1)tt//3 and (3 = l/k B T. The thermo- 
dynamic potential thus takes the form, 

fi = i^TxlnG(k,iw m )+fo, 

P km 

mA^ v [ (2 2 ,i/2 A 2 " 

k 

+ ^E[ ln /(-^+)+ln/(- J Bk-)], (3) 
P k 

where / (x) = [exp((3x) + is the Fermi distribution 
function. We shall confine ourselves to zero temperature, 
where the last term in fl reduces to Sk[-Ek+©( — -Ek+) + 
£ k _e(-£ k _)]. 

The mean field treatment presented above provides a 
simplest unified description for the uniform and spatially 
modulated superfiuids. All these phases have to be de- 
termined using the stationary (saddle point) conditions: 
dfl/dA = 0, dfl/dq — 0, as well as the requirement of 
number conservation n = + = —dQ/dfi. 

We now discuss separately the phase diagram in the 
situations where either the field Sp or the population 
imbalance Sn = —dft/dSp is kept fixed. To this end, 



we trace the evolution of all available mean-field solu- 
tions with increasing the dimensionless coupling constant 
r] = 1/kFCt, where Uf = (37r 2 n) 1 / 3 is the non-interacting 
Fermi wave vector, and seek the one with lowest energy 
(not the thermodynamic potential). To gain a physical 
insight of the competing ground states, we show in Fig. 
1 the landscape of at a selected set of parameters. At 
q = there is a Sarma solution situated between the triv- 
ial normal state at A = and the local BCS minimum 
Ao and corresponding to a maxmium of fl as a function 
of A. On the other hand, for large enough field mis- 
match, a spatially modulated pairing (known as FFLO 
phase) is driven with q • k F ~ 5 p. This forms another 
local minimum in the landscape. Interestingly, a saddle 
point solution necessarily emerges in order to separate 
the local BCS and FFLO minima. 

It is worth noting that not all the solutions are stable. 
As follows we mainly focus on the stability against phase 
separation by the criterion dSn/dS^ > 0, which indicates 
the formation of an inhomogeneous mixed state. An- 
other stability criterion that the superfluid density must 
be positive could also be readily examined 0. 

Fixed chemical potential difference. — We present in 
Fig. 2a the interaction-field phase diagram, constructed 
by finding out the state with the lowest free energy 
F = Cl + fm. The general structure of the phase di- 
agram can be understood by considering the BCS and 
BEC limits first. In the BCS limit with infinitely small 
attraction, rj — > — oo, the kinetic energy dominates and 
the Cooper pair formation is limited to the two Fermi 
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Figure 3: (color online). Comparison of free energies of avail- 
able mean-field solutions at coupling constants as indicated, 
with the free energy of the normal gas Fo being subtracted. 
E c = A/"(0)Ag CS /2 is the condensation energy for a symmet- 
ric Fermi gas, Af(0) = mk,F /(2it 2 h 2 ) is the density of state at 
the Fermi surface and Abcs = 8 exp[n /(2k Fa) — 2]. 



surface. For Sfi < 5fi\ = 1/V2Aq, the ground state re- 
mains the BCS state. For Sfii < Sfi < <5/i2 — 0.754Ao, 
the Fermi surfaces may be translationally deformed, in 
order to increase the overlap for pairing. A FFLO state 
with spatially varying order parameter is thereby more 
preferable. The transition from BCS to FFLO states is 
of first order. Finally, for 5fi > Sfi2, the system trans- 
lates continuously into a normal Fermi liquid phase. As 
an example, for r/ = —1 we show in Fig. 3a the numerical 
comparison of free energies of various competing states. 

The ground state in the BEC limit of n — » +00 is also 
known on physical grounds. Because of the strong attrac- 
tion, all the spin down fermions are likely to pair up with 
atoms in the other state, to form a condensate of tightly 
bounded objects in real space. The distortion of Fermi 
surfaces is prohibited, and then the leftover possibilities 
are the BCS pairing and the Sarma state, as confirmed 
numerically in Fig. 3c. The latter state, in this strong 
coupling limit, is a coherent mixture of condensate and 
a remaining Fermi sea of unpaired atoms. It is energeti- 
cally favorable only for Sfi ~ t as to create an unbound 
fermion, where e& = ti 2 /2ma 2 is the two-body binding en- 
ergy. For sufficiently large mismatch Sfi ~ e +2 2 / 3 eF, the 



condensate disappears and the gas becomes completely 
polarized. Transitions among BCS, Sarma and normal 
phases are continuous. 

The phase diagram in the two limits therefore are en- 
tirely different. Around the BCS-BEC crossover one 
could image a qualitatively change. In particular, the 
spatially varying FFLO and saddle point phases should 
cease to exist with increasing the coupling. We find nu- 
merically (i.e., see Fig. 3b) that for 0.15 < 77 < 0.40 the 
system goes from BCS to the normal state, without expe- 
riencing the FFLO nor the Sarma phase. Our mean-field 
finding is in sharp contrast with a previous proposal in 
Ref. p|, where a direct transition from FFLO to Sarma 
phase are anticipated. This anticipation is another topo- 
logical possibility to connect the two limits. 

In Fig. 2b, by re-expressing Sfi in terms of the non- 
interacting Fermi energy, we compare our results of the 
critical Sfi for superfluid-to-normal transitions with the 
quantum Monte Carlo estimate j^] and the recent ex- 
perimental data on the critical Fermi energy difference 
These differences are calculated assum- 
ing a non-interacting dispersion: (SE F /e F ) c = [(1 + 
(<5n/n) c ) 1/3 - (1 - ((5n/n) c ) 1 / 3 ]/2, where [Sn/n) c is the 
measured critical population imbalance (see, i.e., the Fig. 
5 in Ref. 0]). The mean-field prediction is in good agree- 
ment with the Monte Carlo result, but is about two times 
larger than the measurement. This discrepancy should 
not be taken seriously since the mean-field theory is only 
qualitatively valid. On the other hand, only in the weakly 
coupling BCS regime do the chemical potentials equal 
the Fermi energies. Further, a quantitative comparison 
would require the consideration of the external trap. 

Fixed population imbalance. — In this case, the phase di- 
agram is determined by minimizing E = fl + fin + SfiSn. 
As shown in Fig. 4, now the spatially modulated saddle 
point phase and the FFLO phase are energetically prefer- 
able if they exist. Therefore, on the BCS side, with in- 
creasing imbalance the system goes from the saddle point 
state to the FFLO state, and finally turns into a normal 
gas (Fig. 4a). As the interaction strengths increase, the 
FFLO state disappears and the saddle point phase also 
fades away, while the Sarma state starts to be supportive 
(Fig. 4b). In the strong coupling BEC limit, the Sarma 
state becomes the only solution left. 

The above discussion yields a phase diagram in the 
plane of interaction and polarization Sn/n, as plotted 
in Fig. 5. It is topologically similar to the diagram in 
the rf — Sfi plane, except that the BCS pairing phase has 
now been replaced everywhere by the saddle point phase. 
However, it is important to point out that the saddle 
point phase (shadow regions in the figure), together with 
a sliver of the Sarma state, are intrinsically unstable to- 
wards phase separation, since the slope of the plot of Sfi 
versus Sn for these phases is negative, as illustrated in 
the insets of Fig. 4. This is exactly the precursor for 
a spatially inhomogeneous mixed phase Q. Around the 
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Figure 4: (color online). Comparison of energies of compet- 
ing phases. Insets show the chemical potential difference as a 
function of polarization. The arrow in the inset of (b) indi- 
cates a position, above which the slope of the curve becomes 
positive. 
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Figure 5: (color online). Interaction-porlarization phase di- 
agram. The superfluid-to-normal transition boundary (thick 
line) is to be compared with the experimental data (symbols). 
The shadowed region is unstable against phase separation. 



crossover, our prediction for the appearance of the phase 
separation phase is consistent with the experimental ob- 
servations [i|,|2|. 

In Fig. 5 we compare again the predicted boundary 
for the superfluid-to-normal transition with the experi- 
mental findings of critical polarization (5n/n) c [Tj. The 
agreement seems to be qualitatively good. We note, how- 
ever, that the most intriguing FFLO state is not identi- 
fied experimentally. The window for the FFLO state in 
our phase diagram is sizable, but it may shrink rapidly 



with increasing temperature and an external trap as in 
experiments. 

We conclude by discussing the possible effects of quan- 
tum pair fluctuations beyond mean-field. Three remarks 
are in order concerning the r\ — (Sn/n) phase diagram. 
First, though within mean-field the BCS state is strictly 
confined to the horizontal axis (Sn — 0), the inclusion 
of the pair fluctuations may accommodate a finite pop- 
ulation imbalance. As a result, a narrow window for a 
uniform BCS superfluid opens close to the axis of Sn = 
inside the saddle point phase. Secondly, in our mean-field 
theory the phase boundary for the mixed phase is deter- 
mined indirectly from an instability analysis. It can also 
be fixed following the way in Ref. ffj], i.e., by examining 
the energy of an incoherent mixture of some pure states. 
This alternative method requires considering of pair fluc- 
tuations on the strong coupling BEC side. Finally, so far 
we restrict our analysis to the free space. With a finite 
trap one may instead solve the mean-field Bogoliubov-de 
Gennes equations, or, resort to the local density approx- 
imation [1J|. The latter approach is particularly useful 
in order to take into account the pair fluctuations in the 
presence of traps. Details of these issues on quantum 
fluctuations will be presented elsewhere. 
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